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LIST OF SYMBOLS 


Bi Engine damping, lb-ft-sec/rad 

Lag damper coefficient, lb-ft-sec/rad 
e^ Lag hinge offset, ft 

Hp Shaft horse power 

lbs = l£ + 2e^M^ -Hn^e^ , blade moment of inertia about rotor shaft, slug-ft 2 
Ih Moment of inertia of the rotor hub, slug-ft 2 
1 r = Ih + N lbs 

I eq Equivalent polar moment of inertia of engine/drive train system (referenced to rotor shaft 
speed), slug ft 2 

l£ Blade moment of inertia about lag hinge, slug-ft 2 

K s Rotor shaft torsional stiffness, lb-ft/rad 

Lag damper spring constant, lb-ft/rad 
Mass moment of rotor blade about the lag hinge, slug-ft 
m^ Mass of one blade from lag hinge to blade tip, slug 

N Number of blades 

Qa Total aerodynamic torque of all blades about the rotor shaft, lb-ft 
Q^A Aerodynamic torque of one blade about the rotor shaft, lb-ft 

Total torque on a blade about the lag hinge (aerodynamics and lag damper contribution), lb-ft 
Qe Engine torque (referenced to power turbine shaft speed), lb-ft 
r Aircraft yaw rate, rad/s 

r g Engine to rotor nominal rpm ratio 

v Inflow velocity, ft/s 
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w Aircraft body-axis vertical velocity, ft/s 

\|/ Rotor hub angular displacement, rad 

\|/ 1 Gear box side of the rotor shaft angular displacement, rad 

P Flapping displacement (positive upward), rad 

P Flapping rate, rad/s 

C, Lead-lag displacement (positive against rotor rotation), rad 
? Lead-lag rate, rad/s 

Q Rotor rotational speed, rad/s 

0 O Collective pitch input, deg 

5p Pedals displacement, in. 


IV 



SUMMARY 


This paper presents the results of an analytical study conducted to investigate airframe/engine 
interface dynamics, and the influence of rotor speed variations on the flight dynamics of the heli- 
copter in hover, and to explore the potential benefits of using rotor states as additional feedback sig- 
nals in the flight control system. The analytical investigation required the development of a 
parametric high-order helicopter hover model, which included heave/yaw body motion, the rotor 
speed degree of freedom, rotor blade motion in flapping and lead-lag, inflow dynamics, a drive train 
model with a flexible rotor shaft, and an engine/ipm governor. First, the model was used to gain 
insight into the engine/drive train/rotor system dynamics and to obtain an improved simple formula 
for easy estimation of the dominant first torsional mode, which is important in the dynamic integra- 
tion of the engine and airframe system. Then, a linearized version of the model was used to investi- 
gate the effects of rotor speed variations and rotor state feedback on helicopter flight dynamics. 
Results show that, by including rotor speed variations, the effective vertical damping decreases sig- 
nificantly from that calculated with a constant speed assumption, thereby providing a better correla- 
tion with flight test data. Higher closed-loop bandwidths appear to be more readily achievable with 
rotor state feedback. The results also indicate that both aircraft and rotor flapping responses to gust 
disturbance are significantly attenuated when rotor state feedback is used. 


INTRODUCTION 


Dynamic interface problems involving the engine/fuel control and the rotor/drive train/ airframe 
have been encountered in the ground or flight testing of helicopters for a long time (refs. 1-4). They 
are often manifested in severe torsional oscillations in the helicopter rotor drive shaft. In addition, 
problems occur in the fuel control system (refs. 5 and 6), in multi-engine load sharing, excessive 
vibrations, overtorque, and as undesirable rotor speed variations in certain flight conditions such as 
during autorotational recoveries (ref. 7). These problems not only compromise structural integrity 
but also generate adverse effects on handling qualities of the aircraft. Until recently, these dynamic 
interface problems have received little attention by the airframe/flight control system designers 
because they tend to ignore the rotor speed degree of freedom, thus decoupling, by assumption, the 
effects of the engine/fuel control system on the flight dynamics of the helicopter. Although the 
engine/fuel control specialists have paid considerable attention to these problems, they utilize a 
sophisticated engine dynamic model (refs. 8 and 9)-in conjunction with a rather rudimentary model 
for the helicopter rotor/airframe dynamics for the prediction of the overall system performance 
(refs. 6 and 10). As a result, the prediction of the engine/airframe coupling is often inadequate and 
the problems of dynamic incompatibility can sometimes surface later in the ground or flight tests, 
requiring costly modifications to “fix” the problems. 

These problems are further exacerbated by recent design requirements, especially for military 
helicopters, for greater agility and maneuverability. To meet these requirements, a highly responsive 
engine and a highly responsive rotor speed governor are needed.The required expansion in the 
bandwidth of the fuel control system can introduce severe couplings with the lightly damped tor- 
sional dynamic modes of the rotor and drive train systems, compromising the stability margin. 



Compounding the problem is the trend towards using composite rotor blades with attendant lower 
rotor system inertia. With reduced kinetic energy stored in the system, the rotor becomes more sus- 
ceptible to large variations in its rotational speed during rapid maneuvers. These rotor speed tran- 
sients, especially when combined with concern about overtorquing, can severely increase pilot 
workload and lead to underutilization of the maneuver capability of the aircraft. 

Several studies and flight investigations of advanced engine control systems have been con- 
ducted in the past few years (refs. 1 1-13). These control systems generally incorporate lead compen- 
sation through load anticipation logic to enhance the bandwidth of the engine control system. In 
these investigations comprehensive nonlinear simulation models were usually used, thereby provid- 
ing a detailed simulation of the coupled engine/airframe dynamics and the attendant aircraft handling 
qualities. However, the effects of key design parameters on the system stability and performance 
were not well defined. Further, the benefits of using rotor states as additional feedback variables 
were largely left unexplored. 

A recent parametric optimization study (ref. 14) has indicated that significant improvements in 
helicopter agility and maneuverability can be achieved by continuously varying (“C-V”) the rotor 
speed in high performance maneuvers typically encountered in ground attack and air combat mis- 
sions. The C-V rotor control laws were designed, using a point mass representation of the equations 
of motion for the helicopter. Rotor speed was varied continuously as a function of airspeed, com- 
manded load factor, and other lead signals such as control displacements and rates. The rotor speed 
was permitted to vary in the range of 90 to 120% of the nominal values. With such large variations in 
rotor speed, the effects on flight dynamics and the associated flying-qualities implications need to be 
better understood, if the full potential of maneuverability and agility improvements can be realized. 

The objectives of this analytical study are therefore: (1) to take a new look at the airframe/ 
drivetrain/engine dynamic coupling problem using a simplified parametric high-order model, (2) to 
determine, with the simplified model, the extent to which variations in rotor rotational speed influ- 
ence the flight dynamics of the helicopter, and (3) to conduct an exploratory assessment of the bene- 
fits of using rotor states as additional feedback signals in the control system. The study therefore 
began with the development of a parametric high-order helicopter hover model, which includes both 
body motion, rotor/blade degrees of freedom, air mass dynamics, and engine/rotor speed governor 
dynamics. A linearized version of the parametric model was then used to examine the flight/ 
propulsion system integration and associated flight dynamic behavior and to identify potential ben- 
efits of rotor state feedback. 


Simplied Parametric High-Order Hover Model 

The simplified parametric hover model developed and used in this study is based on refer- 
ences 15-18. The model was developed expressly for the coupled vertical/yaw motion of the 
helicopter in hover; therefore, the other four body-degrees-of-freedom (ref. 17) were not included. 
The equation for the vertical motion of the aircraft is identical to that listed in references 15 and 17, 
but the yaw equation is simplified from reference 17 using a stability derivative representation for 
the aerodynamic terms. The engine and rotor speed governor equations, which are similar to those 
used in reference 17, include a simple first order dynamic representation for the engine and a typical 
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proportional-plus-integral rotor speed governor. The drive train dynamic representation includes a 
lumped equivalent inertia of the engine and transmission, engine damping effect, and a flexible rotor 
shaft. The rotor blade dynamics are represented by the coupled flap-lag rigid blade equations of 
motion for the lag-flap-pitch hinge sequence (ref. 18) with co-located flap and lag hinges similar to a 
UH-60 rotor system. A slight modification to reference 18 was made to account for the spring effect 
of either an elastomeric lag damper or the compressibility effect of the hydraulic fluid of a conven- 
tional lag damper. The rotor speed equation was derived by the Lagrange method, using energy 
equations already provided in reference 18. This equation, along with the lead-lag equation, forms a 
“double pendulum” like expression for the rotor hub and the rotor blades. In this study, primarily 
interest is focused on the collective flapping and lead-lag modes in hover; therefore, interblade 
dynamic couplings are neglected. The equation representing the inflow dynamics is identical to that 
shown in reference 15, which considers only the thrust component of the air mass dynamics. In 
summary, this parametric nonlinear model consists of: 

(a) helicopter vertical and yaw motion; 

(b) rotor rotational degree of freedom; 

(c) rotor blade motions in flapping and lead-lag; 

(d) inflow dynamics; 

(e) drive train with flexible rotor shaft; and 

(f) engine dynamic s/rotor speed governor. 

A detailed description of the parametric model is given in reference 19. A linearized state equa- 
tion of the nonlinear model, as shown in Appendix I, will be used later to examine the effects of 
rotor speed variations on flight dynamics and to investigate the potential benefits of rotor state feed- 
back. For now, however, we will first examine the dominant first torsional mode of the engine/drive 
train/rotor system using the nonlinear parametric model just described. 


Estimation of the Natural Frequency of the First Torsional Mode 

The coupled engine/drive train/helicopter rotor has a lightly damped torsional mode, the natural 
frequency of which is typically much less than the rotational frequency of the main rotor. This mode 
is commonly called the first torsional mode (refs. 1 and 2), which involves primarily the opposing 
motion of the engine and the main rotor blades . Because of its relatively low frequency, it can inter- 
act with the fuel control system and at, a high loop gain, can cause system instability or limit cycles, 
as have been well documented in the literature (refs. 5 and 6). A notch filter is usually used in the 
fuel control system (ref. 20) to attenuate this mode. Therefore, a rapid and accurate estimation of the 
frequency of this mode is needed in the preliminary design of the fuel control law. 

Simple models are available in the literature (refs. 2 and 21) to provide a rapid estimation of the 
natural frequency of the first torsional mode. Many years ago, Sanders (ref. 2) proposed a simple 
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model consisting of two masses with two series springs as shown in figure 1 . One of the two 
“masses” represents the equivalent inertia of the engine plus drive train, and the other represents the 
torsional inertia of the rotor blades; the two springs represent the torsional stiffness of the rotor shaft 
and the centrifugal spring of the rotor blades. Sanders shows a good correlation of the prediction of 
the model with the test data of reference 1. More recently, Odder (ref. 21) proposed a simple model, 
which consists of two masses with only a centrifugal spring, assuming a torsionally very stiff rotor 
shaft. Clearly, this is merely a special case of the Sanders model. 

The simplified parametric model described in this paper, will now be used to provide an estimate 
of the natural frequency of the first torsional mode. By focusing only on the engine, drive train, and 
the rotor system (and thus neglecting other rotor dynamics and aircraft motion), a model which con- 
sists of three masses and two springs, as shown in figure 2, is obtained. The governing equations are: 

I eqVl + K s (Vl-V) + Biig 2 lj/i=r g Q E (1) 


where 


Ir V - N {2xj/e c - e ? M$ 2 } = K s (Vi - v) - Q A 


s- 


f e^Mr 

1 + J2__k 


V 


G J 


.. . . 2 Qc 


% 


c 


X R ~ X H +NI bS 
I^S = 1^ + 2e^M^ + m^ej' 

Qa = n Q^,a 

Q? = Qc.a-c c c-^ 


( 2 ) 

(3) 


Equations (2) and (3) reflect the double-pendulum like representation, as mentioned earlier, for 
the hub-blades system with a spring which includes both centrifugal and lag damper spring effects. 
The latter comes from the spring effect of either elastomeric lag damper or the compressibility effect 
of the hydraulic fluid of the conventional lag damper. 

Linearizing the equations (l)-(3) at the reference equilibrium condition, Vi = V = ^ 0 > 

£(o) = £ 0 = 0, £(o) = £ 0 , and retaining only the spring terms result in a set of six state equations. 
This set of six equations consists of three dynamic modes: the rigid body mode, the first torsional 
mode, and the second torsional mode as shown in figure 3. The rigid-body mode, with zero eigen- 
values, represents the mode in which the engine, rotor hub, and the blades rotate together as a rigid 
body, as can be seen by examining the associated eigenvector. The first torsional mode involves pri- 
marily the opposing motion of the engine and the rotor blades. More hub motion will participate in 
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this mode as the torsional stiffness of the rotor shaft increases or the equivalent inertia of the 
engine/drive train system decreases. The second torsional mode involves primarily the hub and 
blades motion. Again, more engine motion will participate in this third mode as the torsional stiff- 
ness of the rotor shaft increases or the equivalent inertia of the engine/drive train system decreases. It 
can be shown (ref. 19) that the natural frequencies of these two torsional modes can readily be calcu- 
lated using the quartic characteristic equation: 

+ djA,^ + d G = 0 (4) 
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Table 1 shows a comparison of three estimates of the first torsional natural frequency of a mod- 
em medium weight helicopter using references 2 and 21, and equation (4). Available test data for the 
aircraft indicate the natural frequency of this mode to be at 17.2 rad/s. The calculation was per- 
formed both with and without considering the spring effects of the elastometric lag damper used for 
this aircraft. It is seen that inclusion of the spring effect of the elastomeric lag damper is very signifi- 
cant. Because of the relatively lower torsional stiffness of the main rotor shaft for this aircraft, the 
estimate from Ockier’s model (ref. 21) is considerably higher than the test value, as expected. The 
estimate from equation (4) appears to be somewhat better than that from Sanders model (ref. 2), due 
to additional consideration of the hub inertia and inertial coupling associated with the double pendu- 
lum formulation. Additional correlations were performed using equation (4) with test data of another 
modem helicopter and the old test data of reference 1. The calculated results for the first torsional 
mode were also found to be in good agreement with the test values. It should be noted however, that 
when the hub inertia is very small, care must be exercised to properly account for inertia contribu- 
tions from various rotor components in equations (l)-(3) to avoid numerically induced instability. 

Equation 4 provides a good insight into the key parameters influencing the lightly damped first 
torsional mode, which, as described earlier, needs to be considered in the design of the engine fuel 
control system. Three major parameters, i.e, torsional stiffness of the rotor shaft, equivalent engine 
inertia, and rotor speed, were examined. Figure 4 shows the effect of the torsional stiffness of the 
rotor shaft on the natural frequencies of the first and the second torsional modes of the engine/drive 
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train/rotor system of the example helicopter (see table 1). As expected, the first torsional frequency 
approaches the estimated value using Odder’ s model, as the torsional stiffness of the rotor shaft 
increases. It is also shown that the percentage increase in the frequency for the 2nd mode is much 
more drastic than for the first mode. The effect of engine/drive train inertia is depicted in figure 5. 
Increasing the equivalent engine inertia (i.e., engine/drive train inertia referenced to the main rotor 
speed to account for the difference in the engine speed and the main rotor speed) results in decreased 
frequency of the first torsional mode, while the frequency of the second torsional mode remains 
essentially unchanged. A reduction in the amount of 22% in the first torsional frequency will result 
for a 100% increase in the equivalent engine/transmission inertia (e.g., due to a growth in transmis- 
sion). The effect of the variations in rotor rotational speed on the natural frequencies of the two tor- 
sional modes is shown in figure 6. Due to the decreased centrifugal spring effects of the rotor blades 
when the rotor rotational speed is decreased, both the frequencies of the two torsional modes 
decrease, with less percentage reduction for the first mode than the second. While the variations in 
rotor speed affect only slightly the two torsional modes (since the reduced centrifugal spring is only 
part of the spring system consisting also of rotor shaft stiffness, and the lag damper spring effects), 
they affect significantly the aircraft response bandwidth and control effectiveness as will be dis- 
cussed next. 


Some Effects of Rotor Speed Variations on Flight Dynamics 

It has been observed for some time now that the effective vertical damping of the rotorcraft in 
hover appears to be considerably smaller than the predicted value based on constant rotor rotational 
speed (refs. 15, 22-24). Flight test experience has also indicated that the helicopter responses, espe- 
cially in the pitch and roll axes, become considerably more sluggish, when the main rotor speed 
droops. Conversely, the aircraft appears to be more responsive as the rotor speed increases. This 
study shows that variations in rotor speed affect significantly the aircraft response bandwidth and 
control effectiveness. For example, as the rotor speed decays, the effective vertical damping (and 
thus the vertical response bandwidth) is reduced significantly. When the rotor rpm is drooped from 
the nominal operating setting, both the roll damping and roll control effectiveness are substantially 
reduced, resulting in a significant reduction in roll response bandwidth and roll control sensitivity. 

By symmetry, the pitch axis also exhibits similar characteristics in hover. 

Figure 7 shows a comparison of the rate-of-climb responses to a step collective input for the lin- 
ear models with constant rotor speed, varying rotor speed with a rotor- speed governor, and varying 
rpm without a governor. These linear models were derived from the the linear model shown in the 
Appendix I. The associated system eigenvalues are shown in table 2. The rotor speed governor is one 
of proportional plus integral type designed for the example helicopter. Its rpm regulation character- 
istics are shown in figure 8. From figure 7 and table 2, it is clear that the effective vertical damping is 
significantly smaller (about 17%) for the varying rpm, with a governor, than for the constant rotor- 
speed case. A comparison of the vertical acceleration responses is shown in figure 9. Again, the 
response for the varying rpm case (with a governor) is considerably different than the constant rpm 
case, both in the initial transient and thereafter. The transient behavior has been shown previously 
(ref. 15) to be a result of the coupled inflow and flapping dynamics. The flattened response following 
the initial transient, which the varying rpm model predicts, is typically seen in flight test data (see 
fig. 10). The flattened acceleration response is a result of the rotor thrust which is shaped by the 
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slow-down-then-speed-up response characteristics of the rotor speed (fig. 8). Note that the predicted 
response can be expected to improve further when allowing for the feedforward effects from collec- 
tive input to the fuel control unit, a feature that the test aircraft of figure 10 had. This is depicted in 
figure 1 1, in which the calculated normal acceleration responses to a step fuel flow are included. 

This flattened acceleration response to collective input, has not been predicted with a constant rotor- 
speed models (refs. 15 and 23). Although improvements in aerodynamic representation were con- 
sidered in reference 23, the effects due to rotor speed variations were not addressed. 

The effect of rpm variations on roll damping and roll control effectiveness was calculated using 
the simplified rotor equations of reference 16. Some results are shown in figures 12(a) and 12(b) 
respectively for the cases with and without collective retrimmed. In actuality, the situation probably 
resides somewhere between the two extreme cases. For both cases, the control effectiveness is 
approximately proportional to the square of the rpm ratio; however, for roll damping, the two cases 
are different as indicated in the figures. In figure 12(b), data from a comprehensive simulation model 
are also shown. This comprehensive helicopter nonlinear simulation model (ref. 25) includes a single 
main rotor model that considers rigid, hinge-restrained rotor blades with flap, lag, and torsional 
degrees of freedom. A three-state Pitt-Peters model for inflow dynamics is also included in this 
model. The quasi-steady values presented in figure 2(b) were reduced, using a low-frequency model 
reduction procedure commonly called residualization method, from high-order linear models 
generated from reference 25 at various rotor speed of interest (from 89% to 1 1 1% of the nominal 
rotor speed of 27 rad/s). These data from reference 25 correlate well with those generated from the 
simplified model of reference 16. The effect of rotor tpm droops on roll rate response to a step lateral 
stick input is shown in figures 13(a) and 13(b), which correspond to the two cases shown in fig- 
ures 12(a) and 12(b). A first order approximation in the roll rate response was used in generating 
these plots. Note that for a first order approximation, the absolute value of the roll damping, Lp, is 
identical to the bandwidth of the roll rate response to lateral control input. It is seen therefore that the 
effects of rpm droops are twofold: (1) reduced bandwidth and (2) reduced control sensitivity. Both 
effects can contribute to the feel of sluggishness perceived by the pilot. Conversely, increased air- 
craft responsiveness may be perceived by the pilot when the rotor speed increases. 


Rotor State Feedback 

The idea of using rotor state feedback to achieve gust alleviation, helicopter attitude stabilization, 
and vibration reduction, has been around for a long time. In earlier days of helicopter development, 
mechanical devices and schemes, such as delta three (to change blade pitch with flapping displace- 
ment), pitch-cone coupling, and Oehmichen (changing blade pitch with flapping displacement of 
adjacent blade) are examples of using rotor state feedback for gust and vibration alleviation. In 1950, 
Miller (ref. 26) first analyzed the effects of feeding back the flapping displacement and flapping 
velocity to cyclic pitch on the stability and control characteristics of the helicopter. In this work. 
Miller also provided a comparative analysis of several mechanical stabilization devices involving use 
of rotor state feedback such as the Young stabilizer bar (ref. 27), the Kaman servo control, and the 
Hiller control rotor (ref. 28), then widely used by helicopter manufacturers. Later, rotor thrust and 
flap moment feedback concepts (refs. 29 and 30) were also proposed by researchers to lower gust 
response and to reduce dynamic rotor loads. 
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Following the introduction of the electronic flight control system for the helicopter in the late 
1950s and early 1960s, rotor state feedback concepts were further pursued both analytically (refs. 31 
and 32) and experimentally involving flight research (ref. 33). Hall and Bryson showed in their 
analytical work (ref. 31) that neglecting the rotor dynamics in the model used to design a high per- 
formance hover autopilot using linear-quadratic-Gaussian (LQG) methodology can result in unstable 
closed-loop response of the more completely modeled system that includes flapping dynamics. Use 
of rotor state feedback involving flapping displacement and flapping velocity eliminates the prob- 
lem. Rotor state feedback is the essential part for the concept of Individual-Blade-Control advanced 
by Ham (ref. 32) for gust and vibration alleviation and attitude stabilization for helicopters. Flight 
investigation of rotor/airframe state feedback was conducted (ref. 33) to assess the feasibility of 
using rotor tip-path-plane motion to provide additional feedback signals to shape rotor and fuselage 
response. However, in these analytical and experimental studies, primary consideration has been 
focused on the flapping and flapping velocity as the only rotor state variables, in addition to the 
usual fuselage states, serving as feedback signals; use of other rotor state variables, such as blade 
lead-lag and rotor speed information, as additional feedback signals has largely been left unexplored. 
In this study, the effects of rotor state feedback are investigated using the high-order parametric 
model developed in this paper. Both flapping and lead-lag rotor states as well as other rotor states 
such as rotor speed and inflow will be considered as feedback signals to explore their potential 
benefits. Comparisons will be made in terms of closed-loop eigenvalues, frequency and transient 
responses for systems with and without rotor state feedback. 


Effect of Rotor State Feedback on Stability 

It is well known (refs. 31 and 34) that, without rotor state feedback (or high-order dynamic com- 
pensation), instability can develop when attempting to greatly enhance the bandwidth of the closed- 
loop system. The bandwidth of the closed-loop system can be enhanced by increasing the relative 
weighting of the state to control using the simple, full state feedback, linear-quadratic-regulator 
(LQR) control-law design method (ref. 31). In reference 31, the study used an analytical model that 
considers constant rotor speed and includes only flapping dynamics in the rotor subsystem. In what 
follows, the LQR method will be applied to the higher-order model developed in this paper, which 
includes variable rotor speed, flap and lead-lag dynamics, inflow, etc. In this exploratory study, 
perfect sensors (with no sensor dynamics or noise), ideal actuators (with no actuator dynamics, posi- 
tion or rate limits), and zero computational delay are considered. To evaluate the effects of achieving 
high bandwidth regulation without using rotor state feedback, two approaches are employed: (1) use 
a quasi-steady model, which assumes instantaneous rotor tilting, for LQR design and evaluate the 
resulting closed-loop system with the high-order model, and (2) use the high-order model in LQR 
design and evaluate the closed-loop system with rotor state feedback gains dropped. 

Table 3 shows the full 10-state and the quasi-steady 5-state models and their respective system 
eigenvalues. The full 10-state model is that shown in the Appendix I, and the quasi-steady 5-state 
model is reduced from the full 10-state model, assuming instantaneous responses in blade flapping, 
lead-lag motion, and inflow variations. As can be seen from table 3, slight adjustments in the 
yaw/vertical modes and the engine/govemor modes of the quasi-steady model are evident, reflecting 
the residualization effect for rotor and air-mass dynamics. The step responses of the rate-of-climb 
and the rpm to collective input are shown in figure 14 for the full and quasi-steady models. As 
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expected, these responses reflect that, as far as the low frequency region is concerned, the quasi- 
steady model matches well with the full-order model. However, with this quasi-steady model and 
with its attendant LQR gains(see table 4), the closed-loop system can become unstable when evalua- 
tion is made with the full-order model as shown in table 5. The instability occurs when the relative 
state to control weighting p increases, thereby increasing the bandwidth of the closed-loop system. 
From the associated eigenvector, it can readily be identified that the unstable mode is the collective 
flapping mode as shown in table 5 and figure 15. Table 5 lists the migration of all the eigenvalues, 
and figure 15 shows only the root locus of the flapping mode. The full-order model with its attendant 
full-state feedback gains increases the closed-loop bandwidth as the relative state to control weight- 
ing increases, as expected. This is also shown in table 5 and figure 15. However, these stable closed- 
loop system configurations become destabilized and eventually become unstable as the relative 
weighting increases, when the rotor state feedback gains are dropped. This is also shown in fig- 
ure 15. Dropping off the rotor state feedback gains is the second approach, as mentioned earlier, of 
attempting to achieve a higher closed-loop bandwidth without using rotor state feedback. It is inter- 
esting to note that, although the two approaches of attempting to enhance closed-loop bandwidth 
without using rotor state feedback are distinctly different, the results are surprisingly similar. 

It is important to assess the sensitivity of certain gains to the characteristics of the closed-loop 
system, since some of the feedback signals may not be readily accessible. It was found that dropping 
off the inflow gain and the integral-of-the rpm-error gain, individually and in combination, produce 
relatively insignificant changes in the closed-loop eigenvalues. Similarly, variations in lead-lag dis- 
placement and rate gains have very little influence on the closed-loop eigenvalues. However, it 
should be noted that this result is specifically for the use of collective lead-lag displacement and rate 
signals in the stability augmentation for vertical/yaw axes of interest in this paper. When augmenta- 
tion in pitch and roll axes is the objective, cyclic lead-lag displacement and rate signals are essential 
rotor states to feedback in order to preserve system stability and achieve high bandwidth control- 
response and disturbance-rejection (ref. 35). In this study, it was found, as expected, that both the 
collective flapping-displacement gain and flapping-rate gain have a dramatic effect especially on the 
frequency and damping of the flapping mode as shown in figure 16. Availability of flapping rate sig- 
nal is therefore essential for feedback to enhance the damping level of the flapping mode. State esti- 
mation methods, either model-based such as a Kalman filter or non-model-based such as a kinematic 
observer (ref. 36) may be required to provide this vital rate signal for feedback usage. 


Effect of Rotor State Feedback on Transient Responses 

The benefits of including rotor states as additional feedback signals, in addition to the usual 
fuselage states, can be seen clearly by examining the transient responses of the closed-loop system. 
Figure 17 shows the effect of rotor state feedback on vertical velocity transient response to a unit 
initial vertical velocity perturbation. As the relative state to control weighting increases (and thus the 
closed-loop bandwidth increases), the vertical velocity response becomes increasingly oscillatory 
when rotor state feedback is not used. Eventually, the closed-loop system becomes unstable as is 
clear from figure 15. With rotor state feedback, the vertical velocity response remains smooth while 
the response becomes faster and faster with increasing feedback gains. This quickened vertical 
velocity response is accomplished, of course, at the expense of the increased flapping response as 
shown in figure 18. However, with damping added to the closed-loop system by the flapping rate 
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feedback, it is seen that the peak flapping response is significantly smaller for the case with rotor 
state feedback than for the case without rotor state feedback. Also, as expected, the flapping response 
remains smooth as the system bandwidth increases. 

The quickened vertical response has been achieved at the expense of the significantly increased 
responses in collective flapping motion. There are also significant increases in lead-lag motion and 
in the rotor speed variations as shown in figures 19 and 20. These responses are strongly influenced 
by the participation of the dominant collective flapping mode that is destabilized when rotor state 
feedback is not used. With rotor state feedback, both the lead-lag and rotor speed responses remain 
smooth as the level of augmentation is increased, in contrast to the oscillatory response in the case of 
no rotor state feedback. It is interesting to note, from figure 20, that, in contrast to the unaugmented 
case, the action of the augmented system is to allow an initial reduction in the rotor speed to provide 
the energy needed to achieve a quick gain in altitude. Thus, exchange of the kinetic energy of the 
rotor system for the potential energy of the aircraft has taken place. This exchange of energy is initi- 
ated by a rapid collective pulse type input via feedback control. This collective input, as shown in 
figure 21, provides a quick augmentation of thrust and causes the rotor speed to droop, before the 
engine begins to supply the power called for by the rotor speed governor. 

At the high gain level of the augmented system without rotor state feedback, the collective input 
becomes very oscillatory due to the destabilized collective flapping mode as discussed earlier. 


Effect of Rotor State Feedback on Disturbance Ejection 

One of the major benefits of stability augmentation is the significant reduction in the aircraft and 
rotor responses to gust disturbances. Figure 22 shows the effect of rotor state feedback on gust 
response. This figure compares the vertical velocity responses to vertical gust input for open loop 
aircraft (the 10 state model) and the closed-loop systems with and without rotor state feedback. It 
shows that the attenuation in gust response with rotor state feedback increases both in magnitude and 
frequency range as the feedback gains increase. At a frequency of 1 rad/s for example, the attenua- 
tion amounts to slightly more than 10 db (3 times reduction) for the moderate gain level (fig. 22(b)), 
with more attenuation in the low frequency region than in the high frequency region. Although gust 
response attenuation is better in the low frequency region for the closed-loop system without using 
rotor state feedback (due partly to slightly higher gains as shown in table 4) compared with that hav- 
ing rotor state feedback, significant amplification takes place in the higher frequency region. 

Because of the destabilization trend with increased gains for the case without rotor state feedback 
(see fig. 15), the amplification in the higher frequency region, in the neighborhood of the flapping 
frequency, becomes greater as the gain level increases. 

Figure 23 shows the effect of rotor state feedback on the flapping response to vertical gust input. 
The attenuation in the flapping response tends to shift to the higher frequency region as the gain 
level increases. Although there is a slight amplification in the very low frequency region (less than 
0.2 rad/s), overall, considerable reduction in the flapping response to the vertical gust input is 
achieved when rotor state feedback is employed. Without rotor state feedback, the flapping response 
of the closed-loop system again exhibits a large amplification in the neighborhood of the flapping 
frequency, due to the destabilization of the flapping mode as the gain level increases. 
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The frequency responses of the lead-lag motion and rotor speed variations to vertical gust input 
are shown in figures 24 and 25. Although small in absolute terms (less than 0.1 deg lead-lag dis- 
placement per 10 ft/s vertical gust in amplitude), some amplification in lead-lag motion takes place 
for the augmented systems in the low frequency region. Part of the increased lead-lag motion is due 
to the increased rotor speed variations in that low frequency region, as can be seen in figure 25. As 
the feedback gains increase, the augmented system involving rotor state feedback achieves a slight 
attenuation in the higher frequency region. The augmented system without using rotor state feedback 
exhibits a substantial amplification in the high frequency region as the feedback gains increase due 
to the destabilized collective flapping mode. 


CONCLUDING REMARKS 


A parametric high-order hover model was developed for the investigation of flight dynamic 
effects of variations in rotor speed and rotor state feedback. The model included vertical/yaw motion 
of the aircraft, rotor speed degree of freedom, flap-lag dynamics, inflow dynamics, a drive train with 
flexible rotor shaft, and engine/governor dynamics. Using this model, an exploratory study was con- 
ducted which showed that 

o With rotor speed variations considered, the aircraft vertical rate and vertical acceleration 
responses change significantly from those predicted with constant rotor speed. In particular, the 
effective vertical damping decreases significantly from that calculated with an assumption of con- 
stant rotor speed. The result thus appears to provide a better match with flight test data. 

o Both pitch and roll control sensitivity and control response bandwidth decrease significantly with 
droops in rotor speed; conversely, control sensitivity and response bandwidth increases with increase 
in rotor speed. This result is in good agreement with effects observed in flight. 

o The simple formula developed in this paper appears to provide an improved estimate of the natu- 
ral frequency of the first torsional mode of the engine/drive train/rotor system; the simple formula 
also provides a good insight into the key parameters influencing this mode, which is important in 
dynamic integration of the engine and airframe system. 

Using the LQR method, effects of rotor state feedback were investigated. The results showed that: 

o With rotor state feedback, a higher closed-loop bandwidth appears to be more readily achievable 
than without using rotor state feedback; without rotor state feedback, instability occurs when 
attempting to increase the bandwidth of the closed-loop system. 

o Another benefit of the rotor state feedback is the rejection of gust disturbances; both aircraft and 
some rotor responses to gust inputs are substantially attenuated when rotor state feedback is used. 

For the vertical gust examined in the paper, some gust response attenuation is achieved in flapping. 
However, amplifications, though small in absolute terms, can occur in such rotor states as lead-lag 
and rotor speed variations. 
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o For the vertical/yaw axes of the aircraft motion considered in this study, the results indicate that 
collective flapping and flapping rate are the two most effective rotor states for feedback; other rotor 
states, such as rotor speed, inflow, and lead-lag are much less effective. It is essential therefore that 
accurate measurements or estimates of the flapping and flapping rate signals be provided. 
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APPENDIX I: LINEARIZED STATE EQUATION 


The ten-state linear equation used in this study, 


x = Ax 4- Bu 


(AI-1) 


where the state x = col. (8Q, 50, 5w, r, 5HP, (3, ^,(3, v), control u= col. (50 o , 5 p ), 


A 

1.0D+03 * 


Columns 

1 thru 

8 


-0.0110 

0.0000 

0.0000 

0.0076 

0.0010 

0.0000 

0.0000 

0.0000 

-0.0052 

0.0000 

0.0000 

0.0025 

0.0000 

0.0000 

0.0000 

-0.0003 

-1 . 1071 

-0.8304 

0.0000 

0.0000 

-0.0118 

0 . 0000 

o.oooo 

0.0089 

-0.0019 

0.0000 

0.0013 

0.0033 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0144 

0.0000 

0.0070 

0.0000 

Columns 

9 thru 

10 


0.0246 

0.0000 



0.0000 

. 0.0000 



-0.6018 

0.0000 



0.0000 

0.0000 



0.0000 

0.0000 



0.02 64 

0.0000 



-0.8157 

-0.0013 



0.0000 

0 . 0000 



0.0000 

0.0000 



0.0000 

-0.0104 




l.OD+03 * 
-0 . 1477 

0.0004 

0.0000 

0.0000 

0.0402 

0.0000 

0.0000 

0.0004 

0.0000 

0.0000 

-0 . 1072 

0.0000 

0.6926 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.8063 

0.0000 


0.0001 

-0.2024 

0.0006 

-1 . 0825 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 . 0024 

-0.0002 

0.0003 

0 . 0000 

0.0000 

0.0000 

0.0000 

-0.0033 

0.0000 

0.0000 

0.0000 

0.0001 

-0.2270 

-0.0024 

-1.2122 

0.0000 

0.0032 

-0.0239 

0 . 0004 

0.0000 

0 . 0000 

0.0000 

0.0000 

o.oooo 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

-0 . 121,4 

0.0000 


was obtained by linearizing the parametric high-order hover model which was configured to simulate 
a utility helicopter similar to a UH-60 aircraft. For this aircraft, the rotor shaft was assumed to be 
torsionally rigid. 


PRECEDING PAGE BLANK NOT FILMED 


17 



Table 1. Comparison of estimates of 1st torsional natural frequency of the example 
helicopter 


Reference 2 1 

Reference 2 

Equation 4 

Without lag damper spring effect 

With lag 

Without lag 

With lag 

Without lag 


damper 

damper 

damper 

damper 


spring 

spring 

spring 

spring 


effect 

effect 

effect 

effect 

23.23 rad/s 

16.23 

15.77 

16.94 

16.23 


Parameter values 
e£ = 2.875 ft 
I H = 120.89 slug-ft 2 
Ieq = 1234.13 slug-ft 2 
l£ - 803.52 slug-ft 2 
= 48.88 slug-ft 
= 7.59 slug 
N -4 

K s = 399,981.60 Ib-ft/rad 
= 35981.00 lb-ft/rad 
£2 = 30.25 rad/s 
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Table 2. Effects of lag damper characteristics and changing rpm on eigenvalues of the unaugmented example helicopter in hover 


Lag damper characteristics 


ZETA, u.c. = 0.7* 


ZETA, 

u.c. = 0.4 

ZETA, 

u.c. = 0.1 

Modes 

Const, rpm 

Varying rpm 

Const, rpm 

Varying rpm 

Const, rpm 

Varying rpm 

Vertical 

"0.29 

-0.24 

-0.29 

-0.24 

-0.29 

-0.24 

Yaw/(rpm) 

-0.34 

-0.32/-0.58 

-0.34 

-0.32/-0.58 

-0.34 

-0.32/-0.58 

Engine 

-3.33 

-3.33 

-3.33 

-3.33 

-3.33 

-3.33 

Inflow 

-14.38 

-14.60 

-14.42 

-14.61 

-14.44 

-14.59 

Collective flapping 

-9.87 ±j 23.11 

-9.71 ±j 23.01 

-9.82 ±j 23.12 

-9.71 ±j 23.01 

-9.78 ±j 23.12 

-9.65 ± j 22.95 

Collective lead/lag 

-4.99 ±j 5.17 

-5.21/-232.27 

-2.85 ±j 6.59 

-9.5 1/-1 27.29 

-0.72 ±j 7.14 

-18.14 ± j 29.85 


^Damping ratio of uncoupled lead-lag degree of freedom. 


Table 3. Full 10-state and quasi-steady 5-state linear models and their eigenvalues 


Full 10-state mode Quasi-steady 5-state model 

(instantaneous rotor response) 

x = Ax + Bu Xj = AqX| + Bqu 


"x] xj = col.(5Q J 5Q 5w r 5HP) 
_ x 2_' x 2 = COl.(C P C P v) 


A Q - A Rll - A Rl2 A R22 A R21 
B Q=B Ri - A R 12 A R22 B R2 


u = col.(^e o £p) 



Ar 11 

1 

& 

< 


r BR i i 


A = 



,B = 




_ A R 2i 

Ar 22_ 


. Br 2. 



Eigenvalues 


-0.29 ±j 0.01 (vertical/yaw) -0.29/-0.31 (yaw/vertical) 

-0.11 (rpm) -0. 1 1 (rpm) 

-1.27 ±j 0.77 (engine/govemor) -1.39 ± j 0.71 (engine/govemor) 

-14.61 (inflow) 

-9.71 ± j 23.02 (flapping) 

-5.767-23 1.97 (lag) 
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Table 5. Effect of rotor state feedback on closed-loop eigenvalues 


LQR Method: J = (x' p R xx x + u' R uu u)dt 

Jo 


R xx = diag([l.42, 0, 0.25, 142, 0, 0, 0, 0, 0, 0]), R uu = I 2 




P = 

0.001 

P = 

0.005 

__ P = 

0.010 

Feedback 

Unaugmented 

Evaluated with 

Evaluated with 

Evaluated with 

Evaluated with 

Evaluated with 

Evaluated with 

gain design 

eigenvalues 

5-state model 

10-state model 

5-state model 

10-state model 

5 -state model 

10-state model 

Design for 

-0.29 

-0.29 

-0.29 

-0.29 

-0.29 

-0.29 

-0.29 

5-state 

-0.31 

-0.68 ±j 0.77 

-0.66 ±j 0.77 

-0.68 ± j 0.79 

-0.66 ±j 0.78 

-0.68 ±j 0.79 

-0.66 ±j 0.78 

model 

-1.11 

-2.41 

-2.47 

-2.29 

-2.21 

-2.28 

-2.19 


-1.40 ± j 0.71 

-4.98 

-3.08 

-11.38 

-4.89 

-16.12 

-5.38 




-21.06 


-27.80 


-31.46 




-4.83 ± j 23.84 


-0.62 ±j 26.23 


1.49 ±j 27.91 




-5.76 


-5.77 


-5.78 




-232.37 


-233.50 


-234.36 



With rotor 

Without rotor 

With rotor 

Without rotor 

With rotor 

Without rotor 



state feedback 

state feedback 

state feedback 

state feedback 

state feedback 

state feedback 

Design for 

-0.29 ±j 0.01 

-0.29 

-0.29 

-0.29 

-0.29 

-0.29 

-0.29 

10-state 

-1.11 

-0.65 ±j 0.77 

-0.67 ±j 0.78 

-0.65 ±j 0.78 

-0.67 ±j 0.80 

-0.65 ±j 0.78 

-0.67 ±j 0.81 

model 

-1.27 ± j 0.77 

-2.33 

-2.35 

-2.18 

-2.11 

-2.17 

-2.09 


-14.61 

-3.73 

-3.07 

-5.65 

-4.83 

-6.25 

-5.33 


-9.71 ±j 23.02 

-16.16 

-20.79 

-20.20 

-27.31 

-23.17 

-30.87 


-5.76 

-10.33 ±j 

-4.98 ±j 23.82 

-1 1.99 ± j 

-0.78 ±j 26,21 

-13.32 ±j 

1.38 ±j 27.91 


-231.97 

23.87 

-5.76 

26.02 

-5.77 

27.58 

-5.78 



-5.77 

-232.20 

-5.87 

-232.91 

-5.74 

-233.47 



-232.02 


-232.24 


-232.50 
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Figure 1. Sanders model. 



Hub/blades double pendulum 


Figure 2. A simple model consisting of a mass, a flexible rotor shaft, and a double pendulum for 
hub/blades. 
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Rigid-body mode 



• Engine, hub, and blades 
rotate as a rigid body 


First torsional 
mode 




• Primarily engine and 
blades motion; 

• Hub motion participates more 
as shaft torsional stiffness 
increases or equivalent inertia 
of engine/drive system 
decreases. 


Second torsional 
mode 



• Primarily hub 

and blades motion; 

• Engine motion participates 
more as shaft torsional 
stiffness increases or the 
equivalent inertia of 
engine/drive system 
decreases. 


Figure 3. Torsional modes of the simplified engine/drive train/rotor model. 



(Baseline) 

Rotor shaft stiffness ratio 

Figure 4. Effect of rotor shaft stiffness on engine/drive train/rotor torsional modes. 
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Percent of nominal equivalent engine inertia 

Figure 5. Effect of engine/drive train inertia on engine/drive train/rotor torsional modes. 



Percent (rpm) 

Figure 6. Effect of rotor rpm variations on torsional modes. 


25 


18 



(a) 



(b) 


Figure 7 . Comparison of rate-of-climb responses to collective input: a) transient responses to a step 
(1 deg) input, b) frequency responses. 
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Figure 8. Comparison of rpm responses to step (1 deg) collective input. 



Figure 9. Comparison of vertical acceleration responses to step (1 deg) collective input. 
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Figure 10. Typical flight test data showing acceleration response to a step like collective input in 
hover (UH-60 test data). 



Figure 11. Comparison of vertical acceleration responses to step (1 deg) collective and fuel-flow 
inputs. 
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Figure 12a). Effect of rpm variations on roll damping and roll control effectiveness in hover 
(without collective retrimmed), b) effect of rpm variations on roll damping and roll control 
effectiveness in hover (with collective retrimmed). 
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Figure 13. Effect of 1 pm droop on roll rate response in hover: a) without collective retrimmed, 
b) with collective retrimmed. 
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0 2 4 6 8 10 12 


Time (sec) 

Figure 14a). Comparison of rpm responses to step (1 deg) collective input, b) comparison of rate-of- 
climb responses to step (1 deg) collective input. 


31 



Design for quasi-static modei/evaluate 

with high-order model 


— Design for high-orderd model/drop off 
rotor state feedback gains 


Im (rad/sec) 



Figure 15. Eigenvalues of the flapping mode with and without rotor state feedback. 
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Figure 16. Effect of flapping and flapping velocity feedback on the flapping eigenvalues. 


32 



(3 (rad) Vertical velocity 


Open-loop system 
(10-state model) 


Without rotor state 
feedback 


— With rotor state feedback 


Relative state to control 

weighting p = 0.0002 p = 0.001 


p = 0.005 



Figure 17. Effect of rotor state feedback on vertical velocity transient response. 


- Open-loop system Without rotor state With rotor state 

feedback feedback 

Relative state to control 

weighting p = 0.0002 p = 0.001 p = 0.005 



Figure 18. Effect of rotor state feedback on transient response of flapping. 
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Relative state to control 

weighting p = 0.0002 p = 0.001 p = 0.005 
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Figure 19. Effect of rotor state feedback on collective lead-lag transient response. 
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Figure 20. Effect of rotor state feedback on transient response of rotor speed variations. 
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Figure 21 Effect of rotor state feedback on collective input due to feedback control (p = 0.005). 
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Figure 22. Effect of rotor state feedback on gust response: frequency response of vertical velocity to 
vertical gust input. 
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Frequency (rad/sec) 


Figure 25. Effect of rotor state feedback on gust response: frequency response of rotor speed 
variations to vertical gust input. 
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